Spatial Structures and Giant Number Fluctuations in Models of Active Matter 
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The large scale fluctuations of the ordered state in active matter systems are usually characterised 
by studying the "giant number fluctuations" of particles in any finite volume, as compared to the 
expectations from the central limit theorem. However, in ordering systems, the fluctuations in 
density ordering are often captured through their structure functions deviating from Porod law. In 
this paper we study the relationship between giant number fluctuations and structure functions, 
for different models of active matter as well as other non-equilibrium systems. A unified picture 
emerges, with different models falling in four distinct classes depending on the nature of their 
structure functions. For one class, we show that experimentalists may find Porod law violation, by 
measuring subleading corrections to the number fluctuations. 

PACS numbers: 45.70.Qj, 74.40.Gh 
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Active matter, collections of interacting self-propelled 
particles, are found in many different contexts. Exam- 
ples include bird flocks bacterial colonies [H, actin 
filaments propelled by molecular motors Q and vibrated 
granular rods and disks In these, the "activity" 

refers to conscious decision making or internally gener- 
ated cellular thrusts in the biological systems, or impulses 
from a vibrating plate for the granular systems. The com- 
bination of activity and interaction can lead to macro- 
scopic order [7- 12 1. However, these systems are far from 
equilibrium and the usual notions of equilibrium phase 
transitions come under unexpected challenges [lj, [lij ]. 
In particular, macroscopic order, and large scale fluctua- 
tions reminiscent of critical equilibrium systems, coexist. 

Depending on their dynamics and symmetries, differ- 
ent active matter systems exhibit macroscopic polar, ne- 
matic, and/or density order. Polar ordering has been 
demonstrated in point polar particle (PP) models @, 
experiments with granular disks [f|, and continuum the- 
ories lU 14 1 . For polar rods (PR), continuum theories 
rule out macroscopicpolar ordering [lfj , and experiment 
on mobile bacteria [2| and simulations of models of po- 
lar rods are in agreement 11, 1(3]. Apolar rods (AR) or 
active nematics have been studied experimentally in 
hydrodynamic theories \yj\ . and in simulation fiol ] 
and exhibit nematic and density ordering. 

The density fluctuations in the ordered state has been 
characterised by the number fluctuations of = (n 2 ); — 
(n)f of particles in a finite box of linear size I, where 
n is the particle number. In active matter systems, 
°f ~ ( n ) a w ith oc > 1, indicating "giant" number fluc- 
tuations (GNF) in comparison to what is expected from 
the central limit theorem. The exponent a has been used 
to infer the long range correlations in the system. In two 
dimensions, for the PP [1 i, EH, and PR @, EH sys- 
tems, it is now known that a = 1.6, and for AR systems 
[IdlEl a = 2.0. 

Consider now an active matter system relaxing to its 



ordered state from an initial disordered state. As the 
density order grows with time t, there is an increasing 
macroscopic length scale C(t) over which there is en- 
hanced clustering. Information about the spatial struc- 
tures in such a coarsening system can be obtained by 
studying the spatial density-density correlation function 
C(r,t) = (p(0,i)p(r,t)) where p{r,t) is the local density 
at point r. Systems relaxing to an equilibrium state typ- 
ically exhibit clean domain formation (l8| . resulting in a 
linear form of C(r,f) = a - b\r\/£(t) for |r|/£(i) < 1, 
known as the Porod law [h|. On the other hand, many 
systems relaxing to an nonequilibrium steady state vio- 
late Porod law due to a hierarchy of cluster sizes. Ex- 
amples include sliding particles on fluctuating interfaces 
[20l ] and freely cooling granular gases [2lj |. 

In this paper, we ask the following. First, we ask 
whether coarsening active matter systems obey Porod 
law. A few studies have addressed this question - dis- 
crete models of active nematics 
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and a recent numer- 
ical implementation of a hydrodynamic polar model (23j 
have shown non-Porod behaviour. A further systematic 
study is necessary, and in this paper we show that the 
Porod law is indeed violated by all the models that we 
study. 

Second, we ask whether the fluctuations that con- 
tribute to GNF are the same as those that cause Porod 
law to be violated. In particular, we ask whether the 
large distance behaviour of C(r, t) can be deduced by 
knowing a. In general, C(r, t) contains more informa- 
tion than of, as the latter is derived from the former: 



of(t) = l d / d d v[C(r,t) (pf 



(1) 



where (p) is the mean density. If the integrand decays 
to zero over a length scale £ <C I, then of(t) ~ l d ~ (n) 
for large I, or a = 1. Since a > 1 for active matter, the 
upper limit of Eq. ([TJ should contribute to the integral, 
implying that non-trivial correlations extend beyond the 
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scale C(t) > I. Hence, one would expect the behaviour of 
C(r, t) near \r\/C(t) w 1 to contribute to of (i) in Eq. ([I]), 
but we will see below many interesting exceptions to this. 
In this paper, we show that different active matter sys- 
tems as well as other noncquilibrium systems studied in 
other contexts, fall is four distinct classes based on the 
relation between their of (i) and C(r,t). In case of the 
first type, the small \r\/C(t) <C 1 behaviour of C(r,i) has 
no bearing on the exponent a. For the other three types, 
it does, albeit in three distinct ways. 

Type 1: We start with PP systems. We first study nu- 
merically the Vicsek model Q , which we denote as PP(1), 
in two dimensions. All particles move with constant 
speed vq. The positions and velocity orientations 9i of 
particle i at time t + At are given by (t + At) = r j (t) + 
Vi(t)At, and 9i(t + At) = arg [E k e M i9 k(t))} + A£(i,t), 
where the summation over k is restricted to those satis- 
fying |rfc — Yi\ < R, and £ is white noise over the range 
(— 7r,7r]. It is known that the system undergoes a tran- 
sition from an ordered state to a disordered state as the 
noise strength A is increased Q ■ 

We choose parameter values for which the steady state 
is polar ordered with no density bands, and study numer- 
ically the density structures in the coarsening regime — 
a typical snapshot of the density clusters in shown in 
Fig. 0Ja). 111 tne tmie regime studied, C(r,t) has no 
directional anisotropy. Hydrodynamic theory predicts a 
length scale C{t) ~ t 5 / 6 @, [Tj]. Interestingly, we find 
that C(r, t) and the corresponding scaled structure func- 
tion S(k,t)/£ 2 [see Fig. [I^b)] shows a data collapse for 
a completely different coarsening length C(t) ~ t 1 with 
7 = 0.25 ± 0.05. To understand the physical origin of 
this length scale, we studied the two eigenvalues Ai and 
A2 of the inertia tensor of the largest cluster. Both of 
these grow as ~ t - 5 [see Fig. Q^c)], implying that the 
radii of the large clusters grow as i 5 , determining the 
length scale C(t). The S{k,t) [Fig. D^b)] consists of two 
distinct power laws with exponents —1.2 for small kC(t) 
and — 2.6±0.1 for large kC(t); the former has been known 
in hydrodynamic theory [Sj, but we highlight the latter, 
signifying violation of Porod law. In real space, the latter 
implies that C(r, t) has a cusp of the form a — 6|r/£(i)| /31 
with Pi = 0.6 ±0.1 for |r|/£(i) <C 1, and a second power 
law |r| -1 ? with i) = 0.8 ±0.1 for \r\/C(t) > 1. Due to this 
crossover, the GNF exponent a, determined from Eq. ((T|), 
depends only on the exponent 77: 
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a = 2- rj/d. 



(2) 



In the coarsening regime, from a direct measurement 
of of , we find a = 1.6 [see Fig. 0Jd)], consistent with 
Eq. @, and measurements in the steady state Q. 

Interestingly, the structure functions in two other po- 
lar models have the same qualitative behaviour. A mod- 
ified version of the PP(1) model was studied in Ref. Q, 
which we refer to as PP(2) model. The rules of the 
PP(2) model are same as those of the PP(1) model except 
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FIG. 1: (a)-(d): Simulation results for the PP(1) model (sys- 
tem size 1024 x 1024, p = 1.0, v = 0.5, and A = 0.3). (a) 
Snapshot of a part of the system at t — 3200 showing the do- 
main structure, (b) Plot of scaled structure function decays 
as a power law with exponents —2.6 [large k£(t)] and —1.2 
[small k£(t)\. (c) The eigenvalues of the inertia tensor for the 
largest cluster grow as ~ i 0,5 . (d) Variance of number of ~ 
(iV) 1 ' 6 . (e) PP(2) model (parameters same as PP(1)): scaled 
structure function decays as a power law with exponents —1.8 
[large k£.(t)] and —1.2 [small kC(t)]. (f) PR model (system 
size 1024 x 1024, p = 1.0, v = 0.5, and A = 0.2): scaled 
structure function decays as a power law with exponents —1.8 
[large kC(t)] and -1.2 [small kC{t)\. 



that the the new positions depend on the new velocities: 
Ti(t + At) = r.i(t)+Vi(t + At)At. Most macroscopic fea- 
tures remain the same as PP(1), except that the steady 
state configurations have density bands Q. In the time 
regime that we study, C(r, t) is isotropic and bands do 
not form. The structure function is shown in Fig. (He) - 
we find, as in PP(1), C(t) - t 25 and the small |k|£(i) 
behaviour implies 77 = 0.8. However, for large |k|£(t), 
S(\t,t) is a distinct power law with exponent —1.8, im- 
plying a 'divergence' (as opposed to a cusp) of C(r, t) for 
small |r|/£(i) as |r/£(f)|" /32 , with j3 2 = 0.2±0.1 — again 
showing non-Porod behaviour. 

Next we studied a PR model defined in in Ref. [ill ]. 
The time evolution of the 9i's in the PR model 
differs from that in the PP(2) model: 6i(t + 
At) = arg\E k s^[coB(0 k (t)-0i(t))]ew(i0k(t))] + 
A£(i,t), and A 6 (— 7r/2, n/2]. We find that the scaling 
of C(t), as well as the shape of S(k, t) of the PR model 
is similar to the PP(2) model [sec Fig. QJf)]. 

In summary, the models PP(1), PP(2) and PR share 
the following common features: They have a coarsening 



3 



length scale C(t) ~ t 25 . At small |r|/£(i) the correla- 
tion functions violate Porod law either as a cusp or as a 
power law divergence. For large |r|/£(i), they exhibit a 
generic second power law decay with exponent r) = 0.8, 
which determines the GNF exponent a — 1.6. Thus, for 
polar models (type 1), the non-Porod behaviour, and the 
GNF, characterize distinct sources of fluctuation. 

Type 2: We study the discrete AR model introduced 
in Ref. [l(|. The 9i's now evolve as follows. The trace- 
less two dimensional matrix Qjfc = (vjV^) — \$jk (with Vj 
denoting the components of the unit velocity vectors) is 
calculated, where the average is done over the particles 
that are in the disk of radius R, centered about particle 
i. If 9 denotes the direction of the largest eigenvector of 
Q, then 9 l (t + At) = 9 + A£(i,t), with A e (-tt/2, tt/2]. 
The positions n(t + At) = r,(t) ± Vj(i + At) At, with 
the signs being chosen randomly with equal probability. 
The steady state of the AR model is characterized by 
nematically ordered bands, however, in the coarsening 
regime that we study they do not arise. Rather very 
interesting cell-like structures — low density zones with 
high density contours — form [see Fig. [5Ja)], whose radii 
increase with time. The data for C(r, t) and S(k,t)/£ 2 
for different times, collapse when plotted against |r|/£(i) 
or |k|£(t), with C(t) ~ t - 5 [see Fig. GOV) an d (c)j. We 
make an independent estimate of C(t) by counting the 
number of cell-like zones, thus measuring the mean cell 
radius R c (t). We obtain C(t) ~ R c ~ t - 5 [see Fig.^d)]. 
For |r|/£(i) < 1, C(r,t) ~ a - b\r/£(t)\' 31 shows Porod 
law violation with a cusp singularity fi\ = 0.45 ± 0.05 
determined from S ~ fc£ -2 ' 45 [see Fig. Efc)] . Unlike the 
PP models, there is no second power law regime in C(r, t) 
for |r|/£(t) S> 1. The above cusp singularity of C(r,t) 
is similar to another discrete model of active nematics 
in two dimensions (see Fig. 5 of Ref. [H]), and a com- 
pletely different model of particles sliding under gravity 
on a fluctuating interface in one dimension (see Fig. 2 of 
Ref. [1(J). Due to the similar functional form of C(r,t), 
the number fluctuation from Eq. JT]) is: 
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The leading order behaviour af « (A) 2 is seen in our 
data for the AR model [see Fig. Hfe)], as well as for the 
sliding particle system [see Fig. |2(f)]. Thus, the scaling 
af ~ (A) 2 may arise in systems other than active ne- 
matics. Even in an ordinary density phase segregating 
system with density 1/2 consisting of domains of equal 
length £(i), the correlation function is 1 — \ r\/C(t) (satis- 
fying Porod law) for |r| < 2C(t), the number fluctuation 
is exactly af = (A) 2 - 4(A) 3 /£(t) w (A) 2 for large C(t). 

We make two important observations related to 
Eq. ((3]): (i) First, the data in Fig.EJe) and (f), and also in 
the experiment of vibrated granular rods [4| , show a vis- 
ible deviation from the leading af = (A) 2 behavior. We 
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FIG. 2: (a)-(e): Simulation results for the AR model (system 
size 1024 x 1024, p = 0.5, v = 0.3, and A = 0.08). (a) Snap- 
shot of a part of the system at t — 3200 showing the domain 
structure, (b) C(r, t) versus r/C(t) showing data collapse, (c) 
The scaled structure function is a power law with exponent 
— 2.5 at large kC(t). (d) The mean cell radius R c (t) ~ £(£)■ 
(e) Number fluctuations of ~ (A) 2 ' , (f) Sliding particle 
model (10 6 lattice sites, particle density 0.5): of ~ (N) 2 ' . 



claim that the subleading term in Eq. ([3]) may account 
for this deviation. We confirmed that — of/ (A) 2 + \a\ in- 
deed scale as (A) 23 (consistent with /3i = 0.45, d = 2) 
and (A) 5 (consistent with (3i = 0.5, d = 1), respec- 
tively, for the data in Fig. [^e) and (f). More remark- 
ably, we fitted the published experimental data Q in 
this way, and concluded that the experimental system 
has a (3i w 0.5 — that is Porod law is indeed vio- 
lated and the exponent is close to the AR model above. 
We, thus, suggest that Eq. ((3]) opens up a new pos- 
sibility for experimentalists — by measuring the sub- 
leading corrections to GNF, they can indirectly mea- 
sure Porod law violation, (ii) Second, we do not see a 
power law ~ |k| -2 at small k for SCk) [see Fig. [2jc)], as 
suggested by continuum theor y Il7l|. The same is true 
for the sliding particle model [20j, as well as the clean 
phase separating system discussed above — for the latter 
S(k,t) = 4C(t) sin 2 (kC(t)/2)/(kC(t)) 2 . Instead, these 
three examples have a divergence ~ C d as k — > 0. This 
indicates that a 2 - l d S(k ->• 0) ~ Z d /|k| d - l 2d - (A) 2 . 

Type 3: A situation different from type 2 would arise 
if C(r,t) has a power law form ~ |r/£(t)|~'' over small 
through large |r/£(t)|. On one hand there will be non- 
Porod behaviour, and on the other, the same exponent 
contributes to the GNF exponent a and is given by 
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FIG. 3: Number fluctuations: (a) the granular gas (L = 10 5 , 
p — 1.0, ro = 0.5, and S — 0.008) showing GNF with exponent 
a = 1.8. (b) ballistic aggregation (L = 10 J , and p = 1.0) 
shows normal fluctuation with exponent a — 1.0. 



Eq. ([2]), provided r\ < d. Although we are not aware 
of an active matter system exhibiting such behaviour, 
another noncquilibrium system namely a freely cooling 
granular gas in one-dimension plj serves as an example 
of this type. Its structure function shows that the |k| 
space exponent is —0.8 2l|, and hence 77 = 0.2. We re- 



visit this model, and calculate the of (t) in the coarsening 
regime. The result is shown in Fig. [3{a). We find a new 
GNF exponent value a = 1.8, consistent with Eq. ([2]). 

Type 4'- Central limit theorem as in a = 1, may hold 
in an interesting situation. This is when dense clus- 
ters (with masses scaling as £(t) d ) appear in isolated 
locations, leading to temporal "intcrmittency" . In this 
case, C(r,t) = £(t) d S(r) + f(\r/£(t)\). Due to the pres- 
ence of the (5-function, the form of the scaling function 
f(\r/£(t)\) is irrelevant in the calculation of erf from 
Eq. ([I}. Thus, a = 1. Such situations arise in ag- 
gregation models: diffusive or ballistic, wherein parti- 
cles aggregate on contact conserving mass. The density- 
density correlation function for the ballistic system has 
been studied in molecular dynamics and also for an equiv- 
alent lattice model For |r/£(i)| > 0, the scaling 
form of the correlation function starts from a low value, 
rises linearly and then saturates, with increasing |r/£(i)|. 
While Porod law holds, C(r,t) also has a term £(t)6(r). 
We measure erf in simulations of the lattice version of the 
model, and clearly see a = 1 as predicted (see Fig. [3fb)). 
We note that the type 4 is distinct from the other types 
in another respect: (N k ) cx (N), for all integer k > 2, a 
consequence of statistics being dominated by the largest 
cluster. We check that this is true for ballistic aggrega- 
tion. For all other cases (types 1-3) (N k ) cx (N)~ fk , where 
the exponent 7 is specific to a system. 

In summary, we studied well known discrete models 
of active matter and some other nonequilibrium systems, 
to understand similarities and differences in their den- 
sity structures. All the active matter systems that we 
studied are shown to violate Porod law, and the non- 
Porod behaviour is quantified by various new exponents. 
We categorised the relationship between spatial density- 
density correlation function and giant number fluctuation 
into four classes. The first one, formed by polar parti- 



cles and rods, shows two distinct scaling behaviour at 
small and large length scales. These models were shown 
to have a new coarsening length scale £(t) ~ t 25 . The 
non-Porod behaviour does not influence GNF. In the sec- 
ond class, examples being apolar rods or particles sliding 
on fluctuating surfaces, we showed that the subleading 
corrections to the number fluctuations may help experi- 
mentalists to detect Porod law violation. Also the known 
discrete models belonging to this class exhibit GNF for 
a different reason than proposed by continuum theory 
of active nematics. In the third class, the scaled corre- 
lation function exhibits power law divergence at small 
scale. An example for this class is the one-dimensional 
freely cooling granular gas, for which we found a new 
GNF exponent a = 1.8. Finally, aggregation models are 
examples of a fourth distinct class, for which the num- 
ber fluctuations satisfy central limit theorem. We hope 
that this study will encourage experimentalists to probe 
density structures in detail in the future. 
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